Quantum effects on the BKT phase transition of two-dimensional Josephson arrays 
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The phase diagram of two dimensional Josephson arrays is studied by means of the mapping to the 
quantum XY model. The quantum effects onto the thermodynamics of the system can be evaluated 
with quantitative accuracy by a semiclassical method, the pure-quantum self-consistent harmonic 
approximation, and those of dissipation can be included in the same framework by the Caldeira- 
Leggett model. Within this scheme, the critical temperature of the superconductor-to-insulator 
transition, which is a Berezinskh-Kosterlitz-Thouless one, can be calculated in an extremely easy 
way as a function of the quantum coupling and of the dissipation mechanism. Previous quantum 
Monte Carlo results for the same model appear to be rather inaccurate, while the comparison with 
experimental data leads to conclude that the commonly assumed model is not suitable to describe 
in detail the real system. 
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Two-dimensional Josephson junction arrays (JJA)ETu 
have been recently described by means of the two- 
dimensional quantum XY model. These arrays are made 
by superconducting islands, where the charge carriers 
(Cooper pairs) interact with each other through the 
Coulomb interaction and move between nearest-neighbor 
islands through the Josephson tunneling mechanism. 
The following action turns out to describe the systcrnoO 
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it is a variant of the well-known quantum XY model and 
therefore a Berezinskh-Kosterlitz-Thouless (BKT) tran- 
sition is expected to occur at some finite temperature: its 
value is affected by quantum effects, whose importance 
depends on the relative weight of the two terms in the 
action. The first one represents the Coulomb interaction 
between Cooper pairs with charge q = 2e and with the 
capacitance matrix 



Cij — Co 



n 



(z Sij - y^j 



(2) 



where Co and Ci = rj Co are the self- and mutual ca- 
pacitances of the islands, and d runs over the vector dis- 
placements of the z nearest-neighbors. The second term 
describes the Josephson interaction with coupling J be- 
tween nearest-neighbor islands (ij), being <fn — <fij the 
phase difference between the ith and the jth supercon- 
ducting island. ^From a quantum mechanical point of 
view the superconducting phase operators <pi are canon- 
ically conjugated to the Cooper-pair number operators 
hi, \d>i, hi 
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The phase diagram of square and triangular lattices 
of Josephson Junctions was experimentally^ investigated 



and compared witfuthc results of quantum Monte Carlo 
(QMC) simulations!] of the above model. 

In this paper we present an analytical study of the 
phase diagram of JJA based on the effective potential 
approach, op. pure- quantum self- consistent harmonic ap- 
proximation! (PQSCHA), where the effect of quantum 
dissipation is also considered. In mesoscopic systems, 
like JJA, environmental effects can modify the physical 
properties of the isolated system. In order to study the 
open system, an additional term to the action ([!]) is thus 
insertecul't3: 
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where the kernel matrix K(u) — {Kij(u)} is a real 
symmetric matrix and, as a function of u, is even and 
periodic, K(u) = K(—u) = K([3—u), and satisfies 

J du K(u) — 0. It contains the whole information about 
the environmental coupling, in particular it is related to 
the classical damping memory functions Jij(t) by 



l7ij(K 
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where v n — 2-Kn/h(3 are the Matsubara frequencies, 
%j(s) means the Laplace transform of Jij(t), and 



K„ = 



hi) 



du K(u) cos v n u 



(5) 



is the nth Matsubara component of the kernel matrix. 

In Eq. (Q), the Coulomb term is like a kinetic energy 
contribution (with C/q 2 as the mass matrix), while the 
Josephson term plays the role of the potential energy. For 
large capacitive coupling the classical limit is approached 
and the system behaves as a classical XY model, dis- 
playing a BKT phase transition^ at the classical critical 
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temperature. In the opposite limit, the energy cost for 
transferring charges between neighboring islands is high, 
so the charges tend to be localized and phase ordering 
tends to be suppressed at lower temperatures. An im- 
portant role in this system is played by the Coulomb 
interaction range, that can be quantified by means of 
the parameter 77, i.e. the ratio between the mutual- and 
the self-capacitance of the junctions. The connection be- 
tween 77 and the charging interaction range can be found 
observing that the lattes is proportional to the inverse of 
the capacitance matrixa. The Josephson coupling of an 
island with its neighbors, Jz, defines the overall energy 
scale for the system, while the characteristic quantum 
energy scale Ml can be identified considering the bare 
dispersion relation, 
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where = 1 — - J^d cos(fe-d), and choosing fl as the 
maximum frequency 



n 2 



q 2 J 
fr 2 Cn 



□ lattice , 



1 + 8*? (?) 
TT¥ - A lattice. 



One can define a meaningful quantum coupling param- 
eter, which rules the importance of 'quanticity' in the 
system, as the ratio between the two energy scales 



(8) 



The thermodynamic properties of the system are stud- 
ied here by means of the PQSCHA, which_was recently 
extended to treat quantum open systemsEj. By means 
of this approximation scheme the thermodynamics of 
the quantum system (|l|) with dissipation (||) can be re- 
duced to an effective classical problem, where the pure- 
quantum part of the fluctuations is taken into account at 
the self-consistent harmonic level through a temperature- 
dependent renormalization coefficient, while nonlinear 
thermal fluctuations and quantum harmonic excitations 
are fully accounted for. Following the prescription of 



Rcf. 12 the PQSCHA effective potential for our system 
reads 



V e s = J cff ^2 t 1 _ cos ( 
<ii> 
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where we do not consider some additive uniform terms, 
since they do influence neither thermal averages nor the 
critical behavior. The parameter J e g includes the pure- 
quantum corrections; it reads 

Jeff (t, 5,^^) = Je-W*™™ , (10) 

where the renormalization coefficient T>{t, g, 77, jij) repre- 
sents the pure-quantum fluctuations of the relative phase 



between islands, <pi — 4>i-\-d, and is a function of the re- 
duced temperature t = (/3 J) , of the quantum coupling 
parameter g, of the Coulomb interaction range 77, and of 
the damping through jij. In the case of a square lattice 
(the generalization to other kind of lattices is straight- 
forward) the explicit expression of the renormalization 
coefficient is 



V = 



(1 + 877)* 
2N 



Mfe 
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where the tilde means that the square frequencies v 2 , uu^ 
and K Uj k are expressed in units of f2, the characteristic 
frequency scale. They read 
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where Ck and K n k are the Fourier transforms of (g) 
and (Q) respectively. The renormalization coefficient T>(t) 
is very sensitive to the range of the charging interaction. 
Indeed, for a fixed value of the quantum coupling the 
pure-quantum fluctuations of the phase represented by 
T>(t) are strongly enhanced when 77 increases, and they 
saturate when 77 3> 1/nz. This behavior can be explained 
by writing the dispersion relation of the linear excita- 
tions, 



q 2 J c s z\i k 
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as 77 increases, a larger and larger region of the disper- 
sion relation tends to the constant frequency u> B , and 
the relative portion of the Brillouin zone where ujk differ 
significantly from uj e shrinks as (nzr))^ 1 . On the other 
hand, the low-frequency part of the spectrum does not 
contribute significantly to the pure-quantum coefficient 
T>(t) due to the absence of the n = Matsubara term in 
the summation of Eq. ( |TT| ) . 

This behavior of T>(t) also differs from what was found 
in Ref. [l3| (where V is denoted as —2 In go): this results 
from the different choice of the quantum coupling pa- 
rameter. Our choice, Eq. (j^), takes into account the con- 
tributions of the whole capacitance interaction in order 
to determine a meaningful quantum energy scale (0). In 
particular, while the quantum coupling parameter cho- 
sen in Ref. [l3] contains the self-capacitance term only, 
our g varies continuously with 77 and, in the limits of 
small and large 77, it smoothly connects the two quan- 
tum coupling parameters of Ref. |3|. We notice that in 
Rcf. U] the effective Josephson coupling is erroneously 
set to Jeff/ J = (1 + V/2)e- v / 2 - 1 - V 2 /8, while the 
correct low-coupling approximation of Vcffgives Eq. ([!(]) , 
i.e. Jeff/ J ~ 1 — T>/2. Therefore in Ref. [l3| the quantum 
effects turn out to be significantly underestimated. The 
same mistake is made in Refs. n4l 
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Starting from Eqs. (p[-|Tl]), the quantum thermody- 
namic average of any observable 0(n, 0), e.g. a correla- 
tion function or the helicity modulus, can be expressed 
in terms of its W§yl symbol 0(n, <fi) by a classical-like 
formula given byE3, 

(6(n,4>)) = ((0(n,ct> + t)}) es , (14) 

where 

(•)cff = 2 c 1 |d0(.)e-' 3y -'W (15) 

(Z G = J dct>e- fjv ° (t ^) and (£>(n, + denotes the 
average over the two independent Gaussian distributions 
of n and the former accounts for the total fluctua- 
tions of the number variables, while the latter smears 
the function 0(n, <p) on the scale of the pure-quantum 
fluctuations £ of the phases. 

Within the PQSCHA approach, the BKT critical tem- 
perature i BKT of the system, can be-easily evaluated by 
the following self-consistent relatiorJij: 

t BKT (9,V,lii) = ^bkt e-i^WM.7«) t (16) 

where is the BKT critical temperature of the classi- 
cal XY model, t^ T = 0.8950 for the square lattice and 
^bkt = for the triangular one. 
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FIG. 1. Phase diagram for the quantum XY model, 
Eq. (jl]) , on a square lattice in the g 2 — t plane at fixed value of 
r\ — 100. The line is our result for the undamped system, i.e. 
r = Ti = 0. The circles are the QMC results of Ref. || 

Turning to our results, let us consider firstly the be- 
havior of the undamped system (K n ,k — 0). Our phase 
diagram (i BKT vs. the square coupling g 2 ) for the square 
lattice is compared with quantum Monte Carlo simula- 
tion^ in Fig. [jj We observe that our approach starts, by 
construction, from the correct classical value and 

' BKT 

remains valid as long as D < 1: this means that the 
reduction of the transition temperature due to the quan- 
tum effects can be considered reliable if less than about 



30% of tg^ T - On the other hand, the limiting value pf 
the QMC data for vanishing quantum coupling is 0.9430, 
displaying a significant disagreement with the by now es- 
tablished value 0.895 quoted above for the classical XY 
model. 

In Fig. ^ the transition temperature versus g is com- 
pared with the experimental data both for square and 
triangular lattices. It must be noticed that the a pri- 
ori knowledge of the model parameters (J, Co_ and rj) 
describing the experimental setup is rather poorLl: there- 
fore there is an uncertainty both on the vertical (t = T / J) 
and on the horizontal [g as given by Eqs. (|| and 0)] scale. 
Indeed, the values of J used in Ref. fj] give extrapolated 
classical transition temperatures which are larger than 
the known theoretical ones, both for the square (0.95 
against 0.895) and the triangular (1.7 against 1.36) lat- 
tice. In order to avoid such systematic error, each set of 
data in Fig. |^ is normalized to the corresponding classical 
extrapolated value, as already done in Ref. 0. The over- 
all agreement appears to be rather good even though the 
experiments present a more rapid decrease of the tran- 
sition temperature for increasing values of g: this could 
suggest that the pure XY model is insufficient to explain 
this behavior. 
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FIG. 2. Phase diagram for the quantum XY model, Eq. (Q) , 
on triangular (on the left) and square (on the right) lattice 
in the g — t plane at fixed value of r\ = 100. Solid line: 
ro = Ti = 0; dashed line: Fo = 3 and Ti = 0; dot-dashed 
line: To = and Ti = 3. The black squares and triangles are 
the experimental data of Ref. 

Let us now introduce the dissipation. Usually, the 
damping is described in terms of shunt resistors connect- 
ing the islands to ground or/and shunt resistors in par- 
allel to the junctions^. Furthermore the environmental 
coupling is taken to be of Ohmic type, i.e. 

7ij0) = Jij = ^RqGij (17) 

does not depend on s. Rq = 2nh/q 2 is the quantum re- 
sistance and Gij is the conductance matrix, which reads 

Gij = -£-8ij + j^( z s ij - Si >i+ d ) > ( 18 ) 
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where i?o and Ri are resistive shunts to the ground and 
between islands, respectively. 

However, a strictly Ohmic damping leads to unphysical 
results as the logarithrniG,Jiltra-violet divergence of the 
fluctuations of moment all!113, i.e. the number of Cooper 
pairs in each island. The simplest way to consider the 
inertia in the response of the dissipation bath is to use 
the Drude modelEJ, that consists of a real-time memory 
damping with exponential decay of the form 



lit) = 0(t) 



1 e"*A 



(19) 



where 9(t) is the step function and r _1 is the Drude ultra- 
violet cutoff frequency; the Ohmic behavior is recovered 
in the limit r <C t. With this regularization Eq. jl7| ) 
becomes 



7ij( s ) = 



Rq 

2tt 
1 

R~i 
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zdi 



1 
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1 



1 + STl 
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The presence of two distinct characteristic times is consis- 
tent with the choice of two independent damping mecha- 
nism, the on-site and the nearest-neighbor one, related to 
Rq and to Ri, respectively. The two characteristic times, 
ti (I = 0, 1), have to be compared with the characteris- 
tic times of the equivalent circuit, obtained generalizing 
the resistively and capacitively shunted junction modeO 
to a 2D array, i.e. RoCo and R\C\. If Tg is smaller 
than RiCi the response of the baths can be considered 
Ohmic; in the opposite case the behavior of the system is 
no longer resistive and the inertia in the response of the 
dissipation bath must be considered. In our calculation, 
we have therefore assumed ti = RgCg as a representative 
values of the unknown characteristic times of the baths. 

In Fig. H we have also plotted the modification of the 
phase diagram due to damping for realistic values of dis- 
sipation: we use as dissipation parameters the dimen- 
sionless quantities, Ti = Rq/Ri. The comparison with 
the experimental data shows that this kind of dissipation 
is not very relevant for low values of g. Increasing g, but 
remaining in the range where our approach is valid, the 
dissipation appears to affect the phase diagram in the 
opposite way with respect to the tendency of the actual 
experiments. Although our results improve the quantita- 
tive accuracy, these agree with the qualitative behavior 
already found in previous worksEJ. 

The problem of a theoretical explanation of the phase 
diagram of JJA is thus open, and it might be worthy 
to investigate whether the common schematization of 
a Caldeira-Leggett coupling trough the phase variables, 
Eq. (||), is fully justified from a fundamental point of 
view for describing a resistive shunt. Moreover, it is also 
known that dissipation does not have a quenching effect 
onto the fluctuations of all dynamical variables, as the 
simple case of the damped harmonic oscillator shows£3, 
and it might be possible to reasonably modify the mech- 



anism of the environmental interaction in order to repro- 
duce the observed phase diagram. 
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